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A  CONTINUUM  DAMAGE  MODEL 


FOR  VISCOELASTIC  MATERIALS 

by 

Y.  Weitsman 
Texas  A&M  University 

Abstract 

This  paper  presents  a  continuum  damage  model  for  viscoelastic 
materials.  "Damage"  is  expressed  by  two  symmetric,  second-rank  tensors  which 
are  related  to  the  total  areas  of  "active"  and  "passive"  micro-cracks  within  a 
representative  volume  element  of  the  multi-fractured  material. 
Viscoelasticity  is  introduced  through  scalar-valued  internal  state  variables 
that  represent  the  internal  degrees  of  freedom  associated  with  the  motions  of 
long-chain  polymeric  molecules.  The  constitutive  relations  are  established 
from  basic  considerations  of  continuum  mechanics  and  irreversible 
thermodynamics,  with  detailed  expressions  derived  for  the  case  of  initially 
isotropic  materials.  It  is  shown  that  damage  causes  softening  of  the  material 
moduli  as  well  as  changes  in  material  symmetry.  The  special  cases  of  uni¬ 
axial  damage  under  uni-axial  stress  and  the  interaction  of  damage  with 
moisture  diffusion  are  also  considered. 
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1 .  Introduction 

Continuum  damage  modelling,  pioneered  by  Kachanov  (1958),  has  become  a 
most  active  subject  of  research  in  the  recent  decade.  A  comprehensive  review 
of  work  in  this  area  was  given  recently  by  Krajcinovic  (1986). 

Thus  far,  most  of  the  modelling  concerned  the  behavior  materials  like 
rock  and  concrete  which  exhibit  brittle  behavior,  or  metals  -  which  creep  upon 
reaching  a  plastic  range.  Very  little  attention  was  paid  to  the  modelling  of 
damage  in  polymeric  materials  that  creep  viscoelastically .  In  a  recent  work 
by  Schapery  (1980),  viscoelastic  stress-strain  relations  were  modified  to 
include  the  effects  of  damage.  However,  since  those  relations  were  not 
derived  from  a  free  energy  and  the  tensorial  character  of  the  damage  variable 
was  not  identified,  the  above  approach  could  not  specify  the  changes  in 
material  symmetry  which  may  occur  in  damaged  configurations. 

Viscoelastic  response  in  the  presence  of  distributed  micro-flaws  occurs 
in  solid  propellants  and  is  likely  to  be  encountered  in  thermoplastic 
composites.  In  contrast  to  rock  and  metals,  polymeric  materials  absorb 
various  kinds  of  solvents,  which  may  damage  the  polymeric  composites  in  a 
variety  of  ways.  In  the  case  of  water  absorption  by  epoxy-based  composites, 
such  damage  was  noted  by  several  investigators  (e.g.  Ashbee  and  Wyatt  (1969), 
Shirrell  et.-al.  (1979),  Drzal  et.al.  (1985),  and  Jackson  and  Weitsman 
(1985)).  Coupling  between  damage  and  diffusion  was  modelled  by  Weitsman 
(1987d)  assuming  elastic  response. 

In  section  2  of  the  present  work  the  choice  of  internal  state  variables 
to  model  damage  is  motivated  by  micromechanical  considerations.  Consequently 
"damage"  is  represented  by  two  symmetric  second  rank  tensors.  This 
representation  resembles  the  form  employed  by  Murakami  and  Ohno  (1981)  - 
although  their  choice  derived  from  different  physical  and  mathematical 
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considerations . 

In  section  3,  viscoelasticity  is  introduced  by  means  of  scalar  valued 
internal  state  variables  which  correspond  to  the  degrees  of  freedom  afforded 
by  the  motions  of  molecular  chains  within  the  bulk  polymer.  The  thermodynamic 
theory  is  then  applied  when  this  set  of  variables  is  present  in  addition  to 
the  tensor-valued  damage  state  variables. 

In  section  the  special  case  of  an  initialy  isotropic  medium  is 
considered.  Viscoelastic  constitutive  relations  are  formulated,  accounting 
for  both  damage  growth  and  contraction. 

Some  specific  circumstances  are  considered  in  section  5.  These  include 
the  case  of  parallel,  uni-directional  damage  and  stress  and  the  interaction 
between  moisture  flux  and  damage  in  the  absence  of  stresses. 

2.  Basic  Equations.  Damage  and  Internal  State  Variables 

When  a  solid  body  absorbs  solvent  from  the  exterior  the  mass  contained 
within  the  material  volume  of  the  solid  is  not  conserved  and  the  solid-solvent 
mixture  is  a  thermodynamically  open  system.  The  basic  equations  for  this 
circumstance  were  developed  by  several  researchers,  including  Prigogine 
(1968),  Biot  (1973),  Rice  and  Cleary  (1976)  and  Kestin  (1979).  A  version  of 
these  equations  was  given  recently  by  Weitsman  ( 1987a, b) .  Accordingly,  the 
"reduced  entropy  inequality"  is  expressed  as  follows*: 

psS  "  Eij°ij  "  pss*  -  (VT)gi  +  "  p’ifi  ~  8ifi  *  ~  0  (D 

In  eqn.  (1)  $  is  the  Gibbs  free  energy,  pg  -  the  mass  density  of  the 

solid,  -  components  of  the  infinitesimal  strain  tensor,  a.  .  -  components 
J  1 J 


*Eqn.  (1)  follows  an  earlier  expression  [Weitsman,  1987a,  eqn  (17)]  upon 

considering  infinitesimal  deformations  and  p  4>  =  p  tl  -  a  e 

s  s  ij  ij‘ 
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of  Cauchy  stress,  s  -  entropy  density  of  the  solid-solvent  mixture,  T  - 
temperature,  q i  -  components  of  the  heat  flux  vector,  g^^  =  aT/ax^  m  -  solvent 
mass,  f^  -  components  of  the  flux  of  m,  y  and  s  -  the  chemical  potential  and 
entropy  of  the  solvent  in  a  reservoir  in  thermodynamic  equilibrium  with  the 
solid-solvent  mixture. 

To  motivate  the  choice  of  damage  variables,  consider  a  statistically 
representative  volume  element  (RVE)  containing  K  microcracks,  which  may  be 
partially  or  completely  closed.  If  the  material  consists  of  several  phases, 
then  consider  their  equivalent  homogeneous  properties. 

Let  denote  the  open  ("active")  surface  of  the  kth  micro-crack,  and  S*^ 
its  closed  ("passive")  surface,  as  shown  in  Fig.  1.  The  "macro"-level 
stresses  and  strains  are  quantities  averaged  over  the  volume  of  the  RVE. 
Assuming  infinitesimal  deformations,  these  "macro"  quantities  are  given  by  the 
following,  well  known  expressions: 


ij 


J  o.  (x)dV  -  x?  (/  (Vi  +  T°x  )dS 
V  J  Q  1  J  J  1 


and 


K 


+ 


Z 

k=1 


J  (T  (k)x  +  T  (k)x  )dS} 

#  1  J  J  1 

SkUSk 


(2) 


V  =  V  l  eij(^)dV  =  2V  {  |  (uinj°  +  ujni°)dS 


1 


*  £  J  (Vi<k)  *  u  n  (k))dS) 

k='  S.US*  J  J 

k  k 


(3) 


In  eqns.  (2)  and  (3)  SQ  denotes  the  outer  surface  of  the  RVE,  T;  denote 
tractions,  nj^  are  components  of  the  outward  unit  normal,  and  u^^  are 
displacements. 
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If  the  tractions  T.j(k)  are  continuous  across  all  surfaces  Sk  and  Sk*,  and 

T-j0  correspond  to  a  uniform  stress  field  a?.  ,  namely  T.°  =  a?.n.(S  ),  then  it 

‘  J  1  "l J  J  0 

is  well  known  that  a..  =  o.°. 

^  J  ^  J 

Let  Aun  and  au$  denote  normal  and  tangential  components  of  displacement 
discontinuity  at  points  on  some  surface  S.  Then,  by  hypothesis,  on  all 
surfaces  Sk*  Aun^  =  0,  while  au$^  occurs  against  the  resistance  of 
frictional  forces.  It  is  therefore  reasonable  to  expect  that  au$  on  the 
surfaces  Sk*  are  much  smaller  than  au$  on  surfaces  Sk. 

Consequently,  it  is  assumed  that,  instead  of  eqn.  (3),  l.  .  is  given  by 

hi  ■  h  ( { <Vj° +  u/-i°>ds 

0  K 

+  ^  J  (u.n(.k)  +  u.n|k))dS}  (3a) 

k=l  Sk  1  J  J  1 


Although  Sk*  no  longer  appear  explicitly  in  eqn.  (3a),  these  surfaces  of 
discontinuity  still  affect  l . .  through  an  implicit  presence  in  u.-. 

1  J  1 

To  amplify  this  point,  consider  linear  elastic  response.  In  this  case 
o^.  and  e .  •  in  the  RVE  can  be  constructed  by  superposition  of: 

(a)  T.j0  =  o?jnj(S0)  acting  on  the  intact  (undamaged)  RVE,  and 

(b)  T^°  =  0  on  SQ,  T.j(k)  =  -o?j  on  Sk,  while  on  $k  AT^k)  =  0,  Aun^k)  = 

0,  and  Au$(k)  are  related  linearly  to  the  tractions.  (k=l , 2 , . . .K) . 

Since  in  part  (b)  all  tractions  across  SkU  Sk  are  continuous,  obviously 
°ij  =  °ij  <r°r  t^ie  suPerP°sed  solution. 

In  view  of  the  linearity  of  the  problem,  all  displacements  are 
proportional  to  a..0.  Hence  for  part  (a) 

J 


u-a)(x)  S°jka  XjOk° 


(4a) 


while  for  part  (b) 


6 


(b) 


Ui  "  pika(-’  sr  s2’  •••  SK’  si  »  s2  '  - SK  )o 


kl 


(4b) 


In  (4a)  Sijka  denote  the  compliances  of  the  undamaged  material. 
Substitution  of  eqns.  (4)  into  eqn.  (3a)  yields 


e .  . 
1J 


(S“ 


jka  +  2V  k=  1  {  ( Pikit  S1  ’  SK’  Si  >  •••  \  >nj 
Sk 


(k) 


P ik*(?’  S1  ’  V’  ....SK")n[k))dS}c 


ka 


Eqn.  (5)  can  be  expressed  formally  as 

eij  =  Sijka  (ST  - SK;  S1 . SK  )oj 


kl 


(5) 


(6) 


where  S^^  are  the  damage  affected  compliances.  However,  eqn.  (5) 
provides  an  additional  insight  by  implying  that  changes  in  material  symmetry 
are  mostly  due  to  the  configuration  of  "active"  cracks,  while  both  "active" 
and  "passive"  cracks  contribute  to  the  overall  softening.  The  above 
assumption  simplifies  the  constitutive  expressions  developed  in  section  4. 
However,  it  should  be  emphasized  that  the  general  development  of  the 
constitutive  formalism  does  not  depend  on  the  assumptions  that  led  to  the 
approximation  expressed  in  eqn.  (5). 

It  seems  that  a  complete  accounting  of  the  response  of  the  RVE  requires 
the  incorporation  of  2K  parameters  which  represent  the  surfaces  and 
(k=1,...K).  Since  the  shapes  of  those  surfaces  are  generally  unknown,  they 
may  be  approximated  by  some  equivalent  flat  surfaces,  namely  by  vectors  D^)  = 
^(k)^(k)  ancj  ^  (k)  _  £  (k)^(k)  (no  sum  Qn  ^  k=1,  ...K).  However,  since  each 
mierocrack  has  two  equal  and  opposite  surfaces,  the  representation  should  be 
independent  of  the  sense  of  n^kt  Following  Spencer  (1962),  this  is 


accomplished  by  representing  each  microcrack  by  the  symmetric  dyad 
A(k)  =  D(k)D(k).  Similarly  we  have  A*(k)  =  D*(k)D*(k). 

Due  to  the  paucity  of  detailed  information  regarding  the  size  and 
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location  of  the  surfaces  Sk  and  Sk*  it  is  proposed  to  account  for  their  joint 
effect  on  in  eqn.  (4b)  and  in  eqn.  (6)  by  means  of  their  dyadic 
sums,  namely  through 


and 


A  =  I 
*  k  =  1 


* 

A 


I  D*(k)  D*(k) 
k=1 


(7) 


Hence,  the  effects  of  "damage"  are  represented  by  two,  "macro"-level, 

internal  state  variables  that  are  symmetric,  second-rank  tensors. 

The  selection  of  A  to  represent  damage  was  recently  employed  by  Weitsman 

(1987c).  It  was  also  alluded  to  in  an  earlier  work  by  Kachanov  (1980). 

* 

The  tensors  A  and  A  can  be  non-dimensionalized  through  division  by  a 
characteristic  area  A,  e.g.  any  one  of  the  "walls"  of  the  RVE .  This  leads  to 
the  non-dimensional  quantities 


2  *  *  ? 
a  =  A/A  ,  a  =  A  /A 


(8) 


Although  the  present  work  considers  infinitesimal  strains  <<  1),  the 
abovementioned  representation  of  "damage"  does  not  infer  small  a^  and  a^*. 
Consequently,  Taylor  series  expansions  are  useful  in  e  but  not  in  a^  and 


ij' 


As  mentioned  earlier  (Weitsman,  1987c)  the  present  choice  of  damage 
variables  yields,  in  the  case  of  randomly  oriented  micro-cracks  of  equal  area, 
the  intuitively  appealing  result  a.^  =  a2  6 .  .  in  both  two  and  three 
dimensions.  This  outcome,  which  implies  that  randomly  distributed  damage 
affects  the  material's  response  as  a  scalar  quantity  devoid  of  orientation,  is 


not  borne  out  when  a  vectorial  representation  D  =  Dn  is  chosen  for  "damage". 
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V iscoelastic  behavior ,  which  occurs  most  commonly  in  polymeric  materials, 
has  been  attributed  to  the  various  degrees  of  freedom  afforded  by  the  motions 
of  intertwined  long-chain  molecules.  Such  considerations  led  Rouse  (1953), 
Bueche  (1954),  and  Zimm  (1956),  to  express  the  motions  of  molecular  chains 
consisting  of  N  segments  in  forms  which,  upon  transformation  to  principal 
coordinates  and  normal  modes,  resulted  in  N  scalar-valued  internal  state 
variables.  These  formulations  motivated  subsequent  researchers  (Biot  1954a, b, 
Schapery  1964,  1966,  1969)  to  derive  a  theory  of  viscoelasticity  from 
fundamental  concepts  of  irreversible  thermodynamic  processes,  with 
irreversibility  stemming  from  N  internal  state  variables  whose  "growth  laws" 
derived  from  the  motion  of  chain  segments  subjected  to  viscous  resistance. 

It  is  worth  noting  that  in  the  abovementioned  works  the  tensorial 
character  of  the  internal  state  variables  was  not  specified  and  stress-strain 
relations  for  viscoelastic  materials  of  various  symmetries  were  inferred  from 
analogies  with  the  elastic  case.  This  procedure  will  not  carry  over  to  the 
present  work,  which  incorporates  "damage"  by  means  of  two  symmetric  second- 
rank  tensors.  As  shown  in  sections  3  and  4,  the  derivation  of  constitutive 
expressions  and  symmetry  relations  hinges  on  the  identification  of  the 
viscoelastic  internal  state  variables  as  scalars. 

3.  Constitutive  Relations:  Thermodynamic  Considerations. 

Consider  a  viscoelastic  material  with  continuously  distributed  damage 
represen td  by  the  symmetric  second  rank  tensors  a^,  aii*  defined  in  section 
2.  Furthermore  let  (r=l,...N)  denote  N  scalar-valued  internal  state 
variables  which  represent  the  internal  degrees  of  freedom  of  molecular  motion 
within  the  viscoelastic  polymer. 

The  subsequent  formulation  is  based  upon  the  premise  that  the  internal 
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molecular  motions,  represented  by  y^,  occur  on  a  dimensional  scale  which  is 
much  smaller  than  the  micro-damage  expressed  by  a^j  and  a.jj*.  This  suggests 
that  the  viscoelastic  retardation  times  are  unaffected  by  damage,  although 
damage  growth  rates  are  influenced  by  the  inherent  viscoelasticity  of  the 
material . 

io  maintain  tractability ,  and  focus  on  the  effects  of  damage,  we  shall 
assume  a  common  time  dependence  for  all  compliances.  This  assumption  is  valid 
for  many  isotropic  polymers  which  exhibit  a  constant  Poisson’s  ratio. 

Accordingly,  consider  a  Gibbs  free  energy  $>  of  the  form 


0  =  f(aij>  aij>  aij»  m>  T ) 4> ( m ,  T,  Yr)  (rsl,...  N) 


(9) 


Without  loss  of  generality  let  f  >  0.  Also,  let  <j>  be  dimensionless. 

Consider  a  representative  volume  element  subjected  to  fixed  levels  of 
0ij>  m>  anc^  T.  Following  Coleman  and  Gurtin  (1967),  this  circumstance  will 

trigger  an  irreversible  thermodynamic  process  in  the  material,  causing  the 

# 

internal  variables  a^j,  and  y^  to  drift  towards  their  equilibrium  values 

e  *e  p 

a.  ,,  a.  and  y  . 
ij’  ij  Tr 

According  to  the  extremum  principles  of  thermodynamics  (Callen,  I960)  the 
Gibbs  free  energy  attains  a  minimum  at  a^,  ^ij  y^.  Therefore,  at  those 
values  we  have 


and 


8$  n  8$ 

’  sau 


=  0, 


3$ 

* 

3a.  . 
IJ 


=  0 


(10) 


2p? 

3  ®  „  „  3  4>  3  $  # 

2  ‘vS  *  saij  *  —i-  6rrfiaij 


r  "s 


'r  ij 


2  ? 

1  3  $  .  aS 

^  _  6a.  ,6a,  .  + 


8 y  8a .  . 
r  ij 

o 

*  1  3  $  #  # 

6a,,6a,_  +  -  — * - *  6a.  .6a,  >  0 


2  aaij3akit  V  T  3a.. 3a*  ^  2  3a*  ,3a* 

ij  ka  ij  ka 


(ID 
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Eqn.  (11)  holds  at  a^a,  aij6'  Y^6  ? or  aH  arbitrary  variations 


6a .  . ,  6a . .  and  6y  , 
ij’  ij  Tr 


Expanding  $  about  a^,  a^6,  Yp6  employing  eqns.  (9)  and  (10)  and 


assuming  yr  <  Yr  «  1  we  get 


Ho  ,  m,  T,  a  *  a  "e,  y  )  $e  +  \  f (a.  .,  a.®  a  *e,  m,  T)  -f±-  . 


ij'  '  ij’  ij  ’  T'  '  '  2  -'-ij’  “ij’  -ij 


3 Y  3y 

T  S 


G  Q 

(Yr~Yr  )(ys-Ys  )  +  higher  order  terms  in  (y„-y„) 


r  r 


(12) 


Alternately,  an  expansion  about  a.j  =  a^,  a..*  =  a*a  and  Yr  =  0  gives 


*(oij-  »,  T,  a.j,  a.j  ,  rr)  =  Ko1J,  a.j,  a.j  m,  T,  0) 


+  ,  aTl,  a,"  ,  ■.  I,  -S-  rr  .  |  f(V  .jJ,  a  jj,  ..  T)  -fe 


ij »  ij’  ij 

+  higher  order  terms  in  y 


3  Y„ 


y  =0 
r 


r  s 


Y  Y 
r  s 

V° 


(13) 


In  view  of  the  assumption  that  (y^,  y  )  «  1  expansions  (12)  and  (13)  are 
equal,  therefore 


*(°ij,  m,  T,  a.^,  a.  e,  0)  =  H 


+  2  f(oij’  aij’  aij’  m’  T) 

2. 


92  4> 


3 Y  3y 
*r  s 


e  e 

e 


Y  Y 
r  ’  s 


3<t> 


y  -y 

P  P 


3Y, 


3  Y  3  y 
r  s 


V° 


Y  =  Y 
P  P 


and,  in  particular 


(14a) 


(14b) 


3y  3 Y 
r  's 


a'S 


a  3y  9 Y 
Yp=0  r  s 


Y  =  Y 
P  P 


(14c) 


Denote 


3yv 


=  <t>  and  - — f — 
A  r  3y  3y 
Yr=0  'r  's 


V° 


rs 


It  follows  from  eqns.  (9),  (11)  and 


( 1 4a )  that  4>rs  >  0  independently  of  aij  and  ajj*. 
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Following  the  conceptual  models  of  Rouse  (1953)  and  Bueche  (1955),  and 
their  incorporation  into  a  growth  relationship  for  the  internal  state  variable 
as  proposed  by  deGroot  and  Mazur  (1962),  Biot  ( 1954a, b)  and  Schapery  (1964), 
assume  that  ys  depend  only  on  that  portion  of  o  which  concerns  the  internal 
processes  within  the  polymer  (excluding  damage).  Thus  let 


3d)  • 

3 y  =  "  brs  Ys  (r’s  =  (15) 

r 

For  linear  viscoelastic  behavior  brs  =  brs(m,T).  Furthermore,  by 
Onsager's  principle  brs  =  bsp. 

Expanding  <J>(m,T, yr )  in  powers  of  y  ,  truncating  after  the  second  power  in 
view  of  the  foregoing  assumption  that  Yr  <  <  1 ,  we  have 

4>(m,T,Yr)  =  4>(m,T,0)  ♦  4>r  (m,  T,  0)  Yr  +  ~  4*rs(m,T,0)YpYs  (16) 


Substitution  of  (16)  in  (15)  gives 


<t>  + 

r 


)  Y  =  -b  Y 

rs  s  rs  s 


(17) 


Substitution  of  the  specific  form  of  4>,  given  in  eqn.  (9),  into  eqn.  (1) 
yields 


-p 


3f 


3  f 


„  4  o.  .  -  o 

s  3a..  ,  ij  s  3a, 


ij 


ij 


•  3  f 

a .  .  -  p  - j 

ij  Ms  * 


3$ 


3$ 


3a . 


‘ij  Ps  3m  m  “  ps  3T  ^ 


ij 


-P  f 


s  TrV  eij’ij  "  ps  ~  (VT)Si  +  I*  “  -  ~  giV  -  0  (18) 


Since  eqn,  (18)  cannot  be  violated  for  any  process,  familiar  arguments 
(e.g.  Jaunzemis  1967)  give 


£  .  . 
1 J 


•Pj 


and 


3f 

3o .  . 
IJ 


y  =  p. 


3J> 

am 


s  =  - 


3j> 

8T 


(19) 


“p. 


af 

3a.  . 
ij 


a.  . 

ij 


-p. 


3f 


3a. 


•  # 

a.  .-p  f  „ 

IJ  S  3 Y„ 


3  4> 


ij 


V  <qi/T)gi  -  f.^.  -  g.f.  5  >  0  (20) 
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In  particular,  at  fixed  damage,  moisture  and  temperature 
3  4)  •  _ 

’  if  'fr  -  0  (20a) 

r 

Equations  (15)  and  (20a)  give  b  *  *  >  0,  namely  bnc,(m,T)  are 

1 5  r  s  r  o 

components  of  a  symmetric,  semi-positive  definite  matrix.  Since  a  similar 

conclusion  was  reached  earlier  for  <t>  ,  it  is  possible  to  diagnoalize  b  and 

rs  i?  s 

4>  simultaneously  and  rewrite  eqn.  (17)  in  the  form 
r  s 

Br  Yr  +  xr  Yr  =  -<tlr  =  N-  No  sum  on  r)  (21) 

The  solution  of  the  N,  now  uncoupled,  equations  (21)  is 

-t/Z 

Yr  =  K^(1-e  r)  (r  =  1,...  N.  No  sum  on  r)  (22) 


where  =  -4>r/xr  and  Zp  =  Bp/xr  >  0.  Both  Kr  and  Zp  depend  on  m  and  T. 

As  noted  by  Schapery  (1969),  if  the  dependence  of  bn<3  on  m  and  T  is 

I  o 

common ,  namely 

brs(m,  T)  =  b(m,  T)b°s  (b°s  constants)  (23) 

then,  eqn.  (17)  reduces  to 


4>  +  4>  y  +  b°  tt  y  =0 

r  rs  s  rs  dip  Ts 


and,  instead  of  (22),  we  have  the  simpler  result 

-i>/z 


Y  =  K  (1  -  e 
r  r 


) 


(24) 


(25) 


In  equations  (24)  and  (25)  »P  =  t/  b(m,T)  is  the  "reduced  time",  and  b(m, 
T)  is  the  moisture-and-temperature-dependent  shift  factor. 

If,  in  addition,  all  4>rg  =  b^(m,  T ) <p^s  ,  4^  =  constants,  then  eqn. 


(17)  reads 
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dy 

O  O  n 

4>  +  <b  y  +  b  — — 

r  G  rs  Ts  rs  di|» 


=  0 


and  its  solution  reduces  to 


xrbG 


(1  -  e 


-C/Z 


where  5  =  bg(m,  T)  4>  = 


bG(m.T) 

b(m,T) 


(24a) 


(26) 

(27) 


An  expansion  of  <t> (m,T , Yr )  in  powers  of  yr  -  truncating  after  the  second 
power  -  and  employing  the  expression  given  in  eqn.  (26),  yields  the  form* 

Lr  -S/Tr 

*  *  f  bT  O  -  e  >  (28) 

b 

In  the  case  of  a  continuous  spectrum  of  retardation-times  we  obtain 


-t  00 

♦  (».T.t)  -  4-  J  L(t)  (l-e“c/T)d-r  =  4-  D(0  (29) 

Go  dG 

To  proceed  beyond  this  stage  it  is  necessary  to  consider  material 
symmetry.  Isotropy  will  be  selected,  since  it  applies  to  particle  reinforced 
polymeric  composites  as  well  as  provides  the  simplest  illustration. 

4.  Isotropic  Constitutive  Relations 

In  the  case  of  isotropy  f(o^.,  a^.,  a.  .)  depends  on  the  isotropic 
invariants  of  these  three  symmetric,  second-rank  tensors.  Following  Spencer 
(1971)  these  invariants  are: 


ll  °kk’  l2  =  °i  j°ji »  l3  =  °i j°jkaki ’  *4  "  akk*  l5  =  ai jaji *  l6  =  aijajkaki’ 
l7  =  V  l8  =  I5>  *9  =  l6  •  ho  =  aij°ij*  T11  =  aijajk°ki*  !12  =  aij°jkaki* 
!13  =  ai jajk°ks,°8,i *  *14  =  aijaji‘  *15  =  aijajkaki’  *16  =  aijajkaki 


*Note  that  this  form  is  retained  even  for  products  of  Y  because  products  like 
n  -akx  n  -o.x  r 

n  (1-e  )  can  be  expressed  as  sums  i  (1-e  ). 

k=l  k=l 
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A  7% 

*17  =  aijajkakiaai  anc*  n:i-ne  more  invariants  which  include  o.  .  and  a,  (with  or 

1  J  rC  x, 

without  a^p ) .  In  the  above  list  I  =  I(a^j).  The  nine  unlisted  invariants 
are  discarded  in  f  because  it  is  assumed  herein  that,  while  both  "active"  and 
"passive"  microcracks  participate  in  the  softening  of  the  material,  only  the 
active  cracks  affect  its  symmetry.  This  assumption  is  consistent  with  the 

earlier  premises  indicated  in  equations  (3a)  and  (5). 

* 

Expanding  ft  ^ ,  a^,  a^)  in  powers  of  a^j,  terminating  at  the  second 
power,  we  obtain: 


Vkk*  #3aijajk°kl+  5  Vkk»„ 
*  5  ViJ°jl  ’  ft6ai j°ij°kk  *  Vijajk°ki°n 

+  A8aij°ijak!iaam0mk  +  2  A9aij°i jakx°k£, 

2  A 1 0ai ja jkakia&mamn°n2,  +  A11aij°jk°ki 
+  ^i2aijajk  °k  fi,as.  i 


(30) 


In  eqn.  (30)  the  terms  AQ,  A.|,  ...A^  are  functions  of  the  damage 

invariants  1^  -  Ig  and  1^  -  1^,  m,  and  T. 

Strain-stress  relations  are  obtained  from  eqn.  (19)-,.  In  view  of  eqn. 
(30)  we  have 

ea  3psf 

'pq  =  -  —  =  #16pq  +  fi2apq  *  #3aqjajp  *  Vkk5pq  *  Vpq 


pq 


+  A6(aij°ij6pq  +  apq0kk}  +  A7(aqjajp°kk  +  aijajk°ki6pq) 
A8  ^apqdi ja jk°ki  +  ai j°ijaqkakp^ 

A9apqaij°ij  +  A10aq jajpaikaka°ai 
+  A11(aip°qi  +  aqi°iP)  +  A12(aijajp0qi  +  aqjajiV 


(31) 
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The  first  three  terms  on  the  right  side  on  eqn.  (31)  are  related  to 
hygrothermal  expansion  in  the  absence  of  stress.  In  addition,  Alt  A4  and  A5 
should  reduce  to  (cxaT  +  baiti)  and  to  the  "classical"  compliances  Jj_  and  J2  in 
the  absence  of  damage. 


It  can  be  noted  from  eqn.  (31)  that  Furthermore,  if  all 

active  microcracks  are  parallel  to  a  common  plane  it  can  be  shown  that  eqn. 

(31)  corresponds  to  transverse  isotropy  in  that  plane.  Similarly,  if  all 

active  microcracks  are  parallel  to  two  or  three  mutually  perpendicular  planes 

eqn.  (31)  corresponds,  as  expected,  to  orthotropic  symmetry. 

★ 

For  fixed  a^j,  a. o..,m,  and  T  the  viscoelastic  strain-strain  relations 
are  simply  a  product  of  expressions  (29)  and  (31),  namely 


/  ±.  \  _  6  &  • 

epq  1  "  epq 


D(0/bG 


(32) 


For  fixed  damage  but  for  fluctuating  o..,  m,  and  T  the  straightforward 

1  J 

employment  of  superposition  integrals  gives 


>G(n,(t),  T( t) )  epq(t)  =  (alSpq  +  c2apq  ♦  «3aq.ajp)  J'dU-O  d, 
+  <6l‘pq  +  Vpq  +  e3aqjajp>  /Wf)  d. 


pq 

♦  A4«pq  fWs'l^d,  +  a5  fou-v)  d,  *  ... 


So.. 


3a 


A12^ijajp  /  ~ dx  +  aqjaji  l  D^-^)  “I?  dT) 


(33) 


o  "  ^  J '  o 

In  eqn.  (33),  a.  and  f$.  (i=l ,2 ,3) ,  as  well  as  in  all  A.-( j=4,  ...12), 
depend  on  the  invariants  I4  -  I9,  I14  -  I 17  and,  possibly,  on  m  and  T.  Also 

t 


5  =  5(t)  =  J 


du 


b[m(u),  T(u) ] 


and  e'  =  c(t) 


Regarding  the  case  of  fluctuating  damage,  consider  first  the  situation  of 
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monotonic  growth  of  all  micro-cracks.  In  this  circumstance  all  Sk*  =  0  in 
eqn.  (6)  and  the  compliances  are  S.^  =  S.^  (S^t),  S2(t),  ...  SK(t)).  In 
the  context  of  the  continuum  damage  variables  employed  herein  we  have 

S-j jkgr(t)  =  si jka  Under  constant  o.^,  m,  and  T,  and  for  the  elastic 
case,  eqn.  (6)  reads 


which  can  be  expressed  as 

e1j(t)  "  Sijrs‘apq(°»  Crsmn<V,<°»  WWt>)ok?  <*0 

Eqn.  (34),  where  =  S.  ,  shows  that  the  case  of  growing  damge  can  be 

viewed  formally  as  that  of  fixed  damage  with  an  appropriately  modified,  time- 
varying  stress. 

In  view  of  the  assumed  monotonic  growth,  upon  postulating  continuity  of 
Si jka  in  apq »  the  entire  history-dependence  of  Sijk£  is  given  by 

S1jkJt(apq<t>’  ■  Sijkt(ap,(°»  *  l  3SiJaCPq)(T)  dT  (35) 

The  circumstances  considered  in  eqns.  (34)  and  (35)  permit  the  employment 
of  the  correspondence  principle  (Lee,  1955). 

Consequently,  with  a^.  =  a.j(t),  the  strain-stress  relations  take  the 

form 


bG(m,T)  ei j(t)  -  sijka(anin(0))  J  D(5-5')  dr  + 

t  .  t  aS...  (a  )  aa  (u) 

JDU-C)^  lokt<T)  J  - - IS -  duldx  (36) 

o  o  mn 

Upon  identifying  S.jk£  (apq)  with  the  compliances  and  forms  for  the 
intially  isotropic  case  in  eqn.  (31),  the  strain-stress  relations,  with 
monotonically  growing  microcracks,  and  time-dependent  a.  .,  m,  and  T,  read 
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bGepq(t)  =  D*  tiT  0  <“l5pq  +  “2apq  +  “3aqjajp> 

+  m  o  (S1spq  +  e2«pq  +  e3aqjajp)  +  (°kk  0  A4)spq 

+  (“pq  0  A5)  +  [0.  .  0  (A6a,j)l6pq  +  (ckk  o  (A^)) 

+  <"kk  0  <fl7aqjajp»  +  <°kt  0  <A7aijajk»spq 

+  <“ki  0  <Vpqaijajk»  +  <”13  0  <A8aijaqkakp» 

+  <“lj  0  <A9apqaij»  +  <»,1  0  <A10aqjajpa1kakl» 

+  <“qi  0  <Allaip>>  +  <“ip  0  <Allaqi»  +  (“qi  0  (A12aijajp>) 

+  (“ip  0  <fl12aqjaji>»  <37> 

t  r 

where  D*  f  denotes  J  D(c-c')  ¥-  dr  and  (g  o  h)  (t)  denotes 

o  3t 

9(0  lh(amn(0))  +JTf-  ^d.].  (38) 

o  mn 

Obviously,  in  eqn.  (37)  all  &.  and  Aj  (i  =  l,2,3;  j=4,5,...12)  depend 
only  on  the  invariants  I4,  I5  and  Ig  of  a^j. 

In  particular,  if  the  material  is  initially  undamaged  then  (amn(0)) 
=  a,  (amn(0))  =  b,  A4(amn(0))  =  J1  and  A5(amn(0))  =  J2  which  denote  the 
familiar  viscoelastic  expansional  coefficients  and  compliances.  In  this  case 
the  initial  values  of  all  other  terms  -  a2,  <*3,  b2,  63,  Ag,  ...  a12  -  vanish. 

When  the  micro-cracks  experience  both  growth  and  closure  the  current 
configuration,  as  given  by  a-jj(t)  and  a..(t),  may  reflect  distinct  growth  and 
closure  histories.  The  viscoelastic  strain-stress  relation  should  distinguish 
among  such  histories. 

Considering  an  initially  undamaged  material,  the  stress-strain  relations 
under  damage  growth  and  closure  read: 


t  3  a. 

be  =  W0)  1 


+  J  "(c-r) 

0  3T 


ka 


(t)  J 

0 


3t 

t (g)  as. . 


i j  k  a  ^  amn 


a  *) 

.  pq 


3a 


mn 


mn 


3U 


du]di 
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x(c)  as ...  .  (a  a  )  3a 

1Jkt  T  pq  -jga-duldx  (39) 

-  o  3a 

pq 

where  x(g)  and  t(c)  represent  the  time  intervals  for  damage  growth  and 
closure,  respectively.  Those  intervals  may  over-lap. 

Reverting  to  the  initially  isotropic  case,  the  result  is  still  given  by 
eqn.  (37),  except  that  expression  (38)  should  be  modified  to  read 


x(g)  3h(a  a  )  3a 

(g  o  h)  (t)  =  g(x)  [ h( 0)  +  J  - - ES_  _jnn 

3a _  3u 


du 


mn 


-  J 

o 


t(c)  3h(a 


mn , 


ea. 


)  3a 


£9 


3a 


3u 


du] 


(39a) 


pq 


5 .  Some  Special  Sub-Cases, 


(a)  Unidirectional  damage  under  uniaxial  stress. 


Let  all  "active"  and  "passive"  micro-cracks  be  in  planes  normal  to  the 
ection,  whereby  a^  =  a  63i ^3 j  and  a.j  =  a  6^.6^. 

In  this  case  eqn.  (31)  reduces  to  the  transversely  isotropic  relation 


e£ 


£Pq  =  J1°kk6pq  +  J2°pq  +  A(o336pq  +  63p63q°kk) 
+  B  63p63q°33  +  CU3p°3q  +  63qa3p} 


(40) 


In  eqn.  (40)  A,  B,  and  C  represent  combinations  of  the  seven  compliances 
Ag,  ...  A^  of  eqn.  (3D.  Also,  all  five  compliances  in  eqn.  (40)  depend  on  a 
and  a* . 

Consider  the  case  of  uni-axial  stress  =  oQH(t)  and,  in  addition, 
assume  m=0  and  AT=0.  In  this  case  the  strains  are 


el 

Ml 

el 

:33 


el  /  t  v 
e22  =  (J1  +  A)°o 

(J.  +  J„  +  2A  +  B  +  2C)a 
2  o 


(41) 


In  view  of  eqn.  (39a),  the  corresponding  viscoelastic  strains  are 
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(t)/o. 


ll^'/uo  "  e22 


'33 


(t)/a0  =  Jj (0,0)D(t)  +  J  O(t-T)[  3F^a'a  ll|4ffl)dT(g) 

*  0  ' 

3Fi(a.a  )  */T/rU 

3d  * 

(t)/o0  =  |J1(0.0)  +  J2(0.0)]D(t)  +  J  D(t-x)  [— ia,a  ^f!(^))d^(g) 

)  aa*/  /c\\ 

*— T?  dT<c>J 


In  eqns.  (42)  Fj  =  +  A,  F2  =  Jj_  +  J2  +  2A  +  B  +  2C,  and  t(g),  x(c) 

denote  times  of  crack  growth  and  closure,  respectively.  These  expressions 
show  that  Poisson's  ratio,  as  inferred  from  -e11(t)/e33(t) ,  can  serve  as  a 
measure  of  the  effects  of  damage  on  material  properties. 


(b)  Coupling  of  damage  with  moisture  transport. 


Consider  the  flux  of  moisture  fi  in  the  presence  of  damage  but  in  the 
absence  of  stress.  In  addition,  let  the  temperature  remain  constant,  i.e.  aT 
=  0. 


The  flux  fn-  will  depend  on  the  gradient  of  the  chemical  potential  zi  = 

★ 

ay/ax.  as  well  as  on  a^  and  a... 

To  derive  the  flux-gradient  relationship  it  is  necessary  to  generate  all 
the  isotropic  invariants  among  the  two  vectors  and  two  symmetric  second-rank 
tensors  f.j,  z^,  ai j  and  a^.*.  According  to  Spencer  (1971,  page  293)  there  are 
seventeen  such  invariants.*  To  maintain  tractability  we  shall  consider  only 
the  seven  first  terms  from  Spencer's  list.  Accordingly 


invariants  of  the  form  U^(n  )..U.  are  irrelevant  for  the  construction  of 

^  ^  J  J 

flux  -  gradient  relationships. 
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i 


a a,  , 
lk  kj 


+  a  ..  a,  .  )z  . 
Jk  ki  j 


(43) 


invariants  Ijj  -  Ig,  1^  -  I^y  listed  in  section  4. 

According  to  eqn.  (43)  the  flux  of  moisture  is  affected  by  changes  in  the 


material  symmetry  caused  by  both  damage  parameters  and  a*iy  Recall  eqns. 
( 1 9) 2 »  (29)  and  (30).  Accordingly,  in  the  absence  of  stress 


"  =  3m  <AoD(t)/bG> 


(44) 


Since  the  boundary  condition  for  the  moisture  transport  process  is 
p(x|x£boundary)  =  yA  (Ambient  vapor),  eqn.  (44)  shows  that  even  under 


constant  ambient  humidity,  the  boundary  condition  will  be  time-dependent.  A 
similar  conclusion  was  reached  by  Long  and  Richman  (1959),  Frisch  (1964), 
Jackie  and  Frisch  (1985),  and  Weitsman  (1987a)  -  even  in  the  absence  of 
damage . 

Consider  now  the  gradients  zi  =  3y/3x..  In  view  of  eqn.  (44)  z1  will 
include  gradients  of  damage,  through  the  dependence  of  AQ  on  the  invariants  1^ 
“  I91  1 1 4  -  1 1 7 •  Obviously,  z^  will  also  contain  moisture  gradients. 


We  have : 


D '  (  4)  t  _3  ,  V  ^0 


(45) 


b~  3m  b  3a,  „ 
G  kil 


where 
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*  *  a2  AoD(^ 

-  W  a"d  Fm  =  71  <_V - )• 

3m  G 

As  can  be  noted  from  eqn.  (45)  Zj_  also  depends  on  time  through  D(0  and 
its  derivatives. 

Combining  equations  (43)  and  (45)  it  is  observed  that  the  flux  of 
moisture  is  reoriented  into  directions  of  increasing  damage.  Since  damage 
growth  can  be  attributed  to  moisture  [Ashbee  and  Wyatt  (1969),  Drzal  et.  al 
(1985),  Jackson  and  Weitsman  (1985),  and  Weitsman  (1987d)],  the  latter  two 
observations  would  indicate  that  moisture  absorption  and  damage  form  a 
synergistic  mechanism. 

6 .  Concluding  Remarks. 


This  paper  presented  a  continuum-damage  model  for  viscoelastic  materials 
under  several  simplifying  assumptions.  The  distinct  micro-cracks  were 
represented  by  two  continuum-level  internal  state  variables.  These  variables, 
which  were  expressed  by  non-dimensional  symmetric,  second-rank  tensors, 
represented  the  total  areas  of  "active"  and  "passive"  micro-cracks  contained 
within  a  representative  volume  element  of  the  material.  It  was  further 
assumed  that  all  viscoelastic  material  responses  follow  a  common  time- 
dependence,  and  all  detailed  expressions  were  derived  for  materials  that  were 
isotropic  prior  to  the  onset  of  damage. 

The  resulting  expressions  show  that,  in  the  presence  of  damage,  the 
response  of  viscoelastic  materials  contains  two  time-dependent  phenomena.  The 
first  is  attributable  to  the  "inherent"  visco-elastic  behavior  of  the 
undamaged  material,  while  the  second  is  due  to  time-dependent  damage  growth. 
In  addition,  damage  induces  changes  in  the  global  symmetry  of  the  material. 

As  can  be  noted  from  eqns.  (37)  -  (39),  the  coupling  of  damage  and 
viscoelasticity  introduces  very  significant  complexities  in  the  material's 
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response  and  substantial  work  remains  to  be  done  -  both  experimentally  and 
analytically  -  to  attain  a  quantitative  and  practical  understanding  of  the 
phenomenon . 


The  most  important  ingredient  which  remains  missing  in  the  present  work 
concerns  the  growth  relations  for  the  damage  parameters  a.  and  a.*.  To  be 
meaningful,  such  relations  should  be  correlated  with  basic  solutions  for  crack 
growth  and  closure  in  viscoleastic  media.  Several  solutions  for  viscoelastic 
crack  growth  are  already  available  at  the  present  time,  as  reviewed  by 
Kaminski i  (1980)  and  outlined  by  Christensen  (1982).  Additional  research  into 
the  fundamentals  of  viscoelastic  fracture  is  necessary. 
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Fig,  1:  A  sketch  of  a  Representative  Volume  Element 
Containing  K  Microcracks,  with  Open  Surfaces 
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